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Abstract 

The Berry phase has found applications in building topological order para¬ 
meters for certain condensed matter systems. The question whether some 
geometric phase for mixed states can serve the same purpose has been 
raised, and proposals are on the table. We analyze the intricate behaviour 
of Uhlmann’s geometric phase in the Kitaev chain at finite temperature, 
and then argue that it captures quite different physics from that intended. 
We also analyze the behaviour of a geometric phase introduced in the 
context of interferometry. For the Kitaev chain, this phase closely mirrors 
that of the Berry phase, and we argue that it merits further investigation. 
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1 Introduction 


Quantum theory offers a rich supply of mathematical concepts, and, in principle, 
each of these raises the question of its physical interpretation. The geometric 
phase is an interesting example of such a concept as new applications of various 
geometric phases continue to appear. We have discussed this circle of ideas 
before, in rather general terms DEI- Here our intention is to see how they play 
out in a context with specihc demands on the interpretation. 

The Berry phase for pure states has found application in the study of topolog¬ 
ically ordered matter, and recently there has been some discussion p m isi El n] 
about the use of geometric phases for mixed states, and the role they can play, 
when topologically ordered matter is kept at hnite temperature. In this paper 
we analyze the behaviour of two geometric phases for mixed states—Uhlmann’s 
phase [HI El [ID] and the inequivalent m interferometric phase [niiiH] — in the 
Kitaev chain, which is a simple model that can undergo a topological phase tran¬ 
sition. For the Uhlmann phase, we hnd a behaviour that is surprisingly intricate. 
However, we will eventually argue that the behaviour is not really tied to the 
physics we want to capture. In contrast to the Uhlmann phase, the interferomet¬ 
ric geometric phase adopts a topological character and makes a discrete jump 
when the multiplicity of the spectrum of the state changes. In our calculation, 
this translates into a jump when the band gap closes in the Kitaev chain. Thus 
it is potentially a useful topological invariant, and we will argue that it deserves 
further study in the context. 

The paper is organized as follows. In Section we introduce the Kitaev 
chain and its associated Gibbs states, and in Sections and we dehne and 
analyze Uhlmann’s geometric phase for these. Section contains a contrasting 
analysis of the interferometric phase, and Section [^contains a discussion and our 
conclusions. 

2 The Kitaev chain 

A brief introduction to topological phases of matter may begin by reminding the 
reader that phases of matter have been very successfully understood by means of 
the concept of a local order parameter, as in Landau’s theory. The new twist is 
that in certain materials the order parameter may take on a global, topological 
character. Indeed, the issue here is whether some particular geometric phase can 
serve as an order parameter for a topological phase Q 

Let us consider a simple model known as the Kitaev chain m in which such 

^The use of the word ‘phase’ in two different senses is unavoidable, and can be traced back 
to Gibbs. 
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behaviour can be seen. It is a 1-dimensional model built from spinless fermions 
+ al^an = ^nm) at ^ sites with the Hamiltonian 


N 

H [-w{aian+i + al+ifln) - l^alun + M(a„a„+i + 

n=l 


( 1 ) 


Here w is the hopping amplitude, /i is the chemical potential and M is known as 
the induced superconducting gap. We assume that the number of sites is large 
and that periodic boundary conditions are used. 

We can use a Fourier transformation to reexpress the Hamiltonian as 

r2'K Ah 

H = I ( 2 ) 

where the Nambu spinor is d'fc = (a^, and 


Hk = —M sin kay 


F , , 

— + tc cos k 




The eigenvalues of this Hamiltonian are 


(3) 


A± = 


± 



+ w cos k 


+ (Msin kY. 


(4) 


Here, to simplify matters, we will set w = M = 1 and m = /i/2. Moreover, for 
convenience, we will represent the Hamiltonian so that the Bloch vector lies in 
the equatorial plane. Thus, we consider Hamiltonians of the form 


Hk 


At 


-Uk ■ cr. 


(5) 


where a = (era,, Uy, cr^) are the Pauli matrices and 
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Hk = — (m + cos k, sin k, 0), 

Afc 

(6) 

Afc = 2\J{m + cos ky + sin^ k. 

(7) 


What is the physics of this? The model has two energy bands (which is 
why we can treat it as a qubit), with a band gap that closes if m = 1. See 
Fig.lD Kitaev found that in an open chain with a finite number of sites, the 
model behaves very differently depending on whether m < 1 or m > 1. In the 
former case, there are two extra states within the gap, which turn out to be 
exponentially localized on the edges of the chain. In the latter case, these states 
are missing. In this sense, there is a kind of phase transition at m = 1. 
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Figure 1: The spectra for as in Eq. Q and m = 0,1/4,2/4 from left-to-right in 
the top row and m = 3/4,1, 5/4 from left-to-right in the bottom row. The fact that it 
matters whether m is larger or smaller than one is not yet apparent. 


In our model, we can discern a topological order parameter by looking at the 
Bloch vector ([^. As we move along the Brillouin zone, parameterized by k, the 
eigenstates move around the equator of the Bloch sphere for m < 1, whereas 
for m > 1 they are rocking back and forth around the east pole. In the former 
case, we pick up a Berry phase factor e™] in the latter case, the Berry phase is 
zero [151 HSl EZ]- For m = 1, the gap closes and the Berry phase is ill-dehned. 
Thus, we can say that there is a topological Z 2 invariant taking different values 
in the two phases, and this invariant serves as the order parameter of the model. 
A systematic classihcation of non-interacting models of this kind is available 

[151 ESI [201 ED- 

After this (absurdly brief) introduction to topological phases of matter, we 
can go on to ask about mixed states. Assume that the chain is in contact with 
a thermal bath at temperature T, and is described by the Gibbs state 

g-r/fe/T I / At \ 

pik) = —-= - 1 -|- tanh -^nk ■ • (8) 

This state dehnes a curve within the equatorial plane of the Bloch ball, as shown 
in Fig. I^for various values of m and T. These curves can be described analyt¬ 
ically. Using polar coordinates on the equatorial plane we have for the angular 
coordinate that 


The case m = 0, which is called the flat band case, is particularly simple because 
the Brillouin zone is mapped to circles centred at the maximally mixed state. 
For m < 1 the curves enclose, for m = 1 the curves pass through and for m > 1 
the curves are conhned to one side of the maximally mixed state. 



ip{k) = - + arctan 


1 — m k 

-tan — 

1 -|- m 2 
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Figure 2: Curves in the Bloch ball for some values of m. In each panel, we see 
the curves swept out by the Gibbs states as we move along the Brillouin zone for 
T = 1/5, 2/5,3/5,4/5,1. The east pole is marked with a dot. 


The question now is if we can define a phase factor for mixed states, analogous 
to the Berry phase for pure states, which can serve as a topological invariant 
describing the phase structure of the Kitaev chain at finite temperature. The 
work by Viyuela et al. iE] at hrst sight suggests that the Uhlmann geometric 
phase |S] can play this role. It would further suggest that there is a sharp 
transition between different phases of the model also at finite T, and that this 
transition happens at an m = mc{T) < 1. For curves in the equatorial plane 
of the Bloch ball the Uhlmann phase factor is indeed a number ±1 which can 
be calculated, and by fixing the curves in this manner Viyuela et al. have added 
an intriguing amount of concreteness to its study. We will add some more, 
and discuss its interpretation afterwards. We will also broaden the perspective 
slightly and study the behaviour of another candidate geometric phase for mixed 
states [12]. A paper with a similar aim has appeared recently [7]. 

3 The Uhlmann phase factor 

The basic idea in Uhlmann’s theory is to let the pure states in the bipartite 
Hilbert space form the total space of a fiber bundle over the mixed states 

on H, see Ref. m. The projection onto the mixed states is V’ '0'^^ = Py 
and the structure group is U{N) acting from the right, "0 —>■ 'iIjU. A geometric 
phase can be associated to any curve in the base manifold once we have defined 
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a parallelism condition for cnrves in the total space—a lift ‘ip{k) of p{k) is said 
to be parallel if for every infinitesimal 6k the probability for the transition from 
to '0(/c + 6k) eqnals the hdelity of p{k) and p{k + 6k), 


^r{^jJ{ky^jJ{k + 6k)) 


Tr p{k)^/‘^p{k + 6k)p{k)^l‘^. 


( 10 ) 


The parallelism condition can also be described in terms of a connection [22] , 
which we here denote by A. In general, no explicit formula for A is known. 
However, Hiibner [23] has derived a formula giving its values along the velocity 
helds of ’square root lifts’ -0 = y/p. The formula of Hiibner reads 


H \ui) 




Pi + Pj 




‘3\i 


( 11 ) 


where the Pi and the \ui) are the eigenvalues and eigenstates of p, and the over¬ 
dot means differentiation with respect to k. It should be pointed out that the 


derivation of (11) assumes that p has full rank, so it fails when the system is in a 


pure state. Indeed, the Uhlmann bundle is dehned for mixed states of full rank 
only. Apart from that, it is insensitive to the spectra of the mixed states. 

Let us postpone a discussion of the physical import of Uhlmann’s construc¬ 
tion and turn directly to an evaluation of the Uhlmann phase for curves in the 
equatorial plane of the Bloch ball. In general, Hiibner’s formula ( [IT| is hard to 
handle, but for the qubit state space it simplihes to 


= (v'ft - i/P2)“(I“i)(“i|“2)(m 2| + |t<2){U2|t<l){t<l|). 


( 12 ) 


Since p depends on k, so do its eigenvectors \ui) and its eigenvalues pi. We are 
interested in curves that stay in the equatorial plane of the Bloch ball, and set 


1 / 1 {pi-p2)e^A 

^ 2 \{pi - P2)e^'^ 1 ) 

and 

l“i) = ^ , l«2) = ^ . 

Then the connection becomes abelian, 

AW ^[~0 ?). 


(13) 


(14) 


(15) 


and no path ordering is necessary in order to compute the symmetry group 
element 


U ik) = exp 


(6kA{'ilj) 


0 \ 

0 e-^^j ’ 


( 16 ) 
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where 


(17) 


A 


1 

2 



dk(p(y^ 


A parallel lift of p(k), then, is 




^ll(k) = ^(k)U(k) 


W (v^+v^)e'^ (v^- 

2 [(^ - (V?h + J • 


(18) 


The unitary l/(k) depends on the special lift we started out with, -0 = y/p, but 
'4’\\{k) does not. Let us remark—since it will be important to us later on—that 
i}\\{k) is in general not periodic, even if the original curve p{k) is closed with 
periodicity 2tt. The function A{k) need not be periodic even if p{k) is. 

Uhlmann’s geometric phase is the argument of the phase factor of the function 

Tr (V^ii (0) V|| (fc)) = ^ (\/pi(0) + V&^) (Vpi(^) + VSW) cos A + 

+ ^ (Vpi( 0) - t/psW) (\/pi(fc) - \/p2W) cos ((p + 71). 

Because we consider curves in the equatorial plane of the Bloch ball, this is a 
real-valued quantity. For small values of k therefore the phase factor equals -|-1, 
and it stays that way until we reach a node on the curve, that is, a point where 
the trace vanishes. At a non-degenerate node the phase factor jumps to —1. If 
the curve was deformed slightly away from the equatorial plane we would see the 
phase building up continuously around the position of the node. 



4 Concrete behaviour of the Uhlmann phase factor 

We now specialize to the case of Gibbs states of the Kitaev chain, for which 


A 


Ai 


Pi = 


(l + tanh^), p 2 = ^(l-tanh^). 


2T' 


( 20 ) 


To simplify matters we consider one of two cases, either the flat band case (m = 0, 
Afc = 2) or the case of closed curves (using Aq = A 27 r)- In either case Pi{k) = 
Pi{0) = Pi, and Eq. (19) simplihes. We hud that 

(v^+v/^)^ = 1+ sech^, iVm-= l-sech^. (21) 


Before inserting this result into Eq. (19), we dehne 

Afc 


X = sech 


2T' 


( 22 ) 
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Recalling that Uhlmann’s phase factor is multiplied by a factor —1 each time 
the curve encounters a node, we see that the task is to hnd those values of k for 
which 


Tr 


) = ^ ((1 + x) cos y4 + (1 — x) COS {ip + R)) = 0. 


In 


this equation, x = x{k;T) as x depends on A*,, and moreover 

A = A{k; T) = - [ dk <p(l — x). 

2 Jo 


(23) 


(24) 


In general this integral is too hard for us to do analytically. 

We hrst consider the simplest case, the flat band case in which = 2, 
X = sech(l/T) and 

(25) 


A = ^(l 
2 ^ 


X) 


The curves are then circles, and the temperature alone determines the purity of 
the states. The limit T —)• 0 can now be taken. When T = 0 the hyperbolic 
secant vanishes, the state is pure, A = ip/2 and the nodes are determined by the 
equation 

ip f o 1' 


cos —h cos — = U 
2 2 




cos 2 


cos^ —-) =0. 

2 2 / 


(26) 


Hence there are nodes at (p = 7r/2, tt, 37r/2. This should be compared with the 
the Berry phase, which has a single node at = tt. However, for a closed circle 
the contributions from the two extra nodes will cancel each other, and we obtain 
the same value for the Uhlmann and Berry phases (as noted by Uhlmann [S]). 

To see what happens as we increase T, that is to say as we shrink the radii 
of the circles, refer to Fig. The hrst thing that happens is that the two 


’extra’ nodes between 7r/2 and vr disappear. Then there is a value of T (called 
the ’critical temperature’ by Viyuela et al. 0) for which the circle coincides with 
the dashed circle. For higher values of T no nodes are encountered during the 
hrst revolution around the Brillouin zone, and consequently Uhlmann’s phase 
factor remains equal to +1. But the evolution of the purihed state 'i/>\\{k) is 
not periodic in k, so a node will be encountered during the second revolution 
(unless the temperature is so high that we are inside the smaller dashed circle), 
see Fig. and so on. 

Fig. 1^ gives the location of the nodes for as many as 5 revolutions. The curves 
dehning the nodes rehect the rather complicated evolution of ipwik). However, 
by inspection we see that some periodicities appear if we consider closed curves. 
The condition for a node to occur after exactly ni revolutions is 


cos (24(27rni)) = cos((l — x)?7,i7r) = 0. 


(27) 





(b) Between one and two turns. 


Figure 3: Where the nodes he if the paths start on the positive x-axis and go counter¬ 
clockwise along circles of constant purity. 


TT 2tt Ztt Att Sn Btt 7 tt 8 9n lOn 



Figure 4: Solutions of Eq. (23) for = 2 and winding numbers up to 5, with x = 
sech(l/T) as the ordinate. For large T, that is for high entanglement, the circles 
traversed in the Bloch ball are small and several revolutions are needed to see a node. 


There are ni solutions for rii revolutions, namely 


X 


ni,n2 


2(ni - ^ 2 ) - 1 
2ni 


ns e {0,1,... ,ui - 1}. 


The node at x = 1/2 repeats after two more revolutions since 

1 

2^1,0 — 3^3,1 — ^5,2 ~ ~ 2' 


(28) 


(29) 


For open curves no special simplihcations are visible. 

Going beyond the flat band case is a more involved story. We hrst ask for 
what value of T we have a node occurring after exactly one turn. This means 
that we have to solve the equation 

1 

cos(A(T))=0, A{T) = - x{k,T)). (30) 

2 Jo 
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(a) The location of the ’critical 
temperature’ Ter = after one 
turn, in agreement with [3]. 


T 


Ts.O 




•h'”.. ^3,1 
Ti^ 


0.0 02 0.4 0.6 08 1.0 1.2 1.4 
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(b) Location of critical tempera¬ 
tures T'rti,n 2 for ni = 1, 2 and 3 
turns. 


Figure 5: Location of critical temperatures T„i,n 2 for ui = 1, 2 and 3 turns. Compare 
with Fig. 1^ which shows the flat band case only. The dashed line in panel (a) gives 
the size of the band gap (with Boltzmann’s constant set to 1). 


This is a straightforward Mathematica calculation, and the result is reported 
in Fig. 5a The resulting curve is in fact a cross-section of any one out of the 


three two dimensional figures given by Viyuela et al. [3] (their Figure 2, which 
is considerably more beautiful than ours). They interpret the resulting values 
of T—less than half of the bandgap at m = 0—as an m-dependent critical 
temperature at which something dramatic happens to the Kitaev chain. We can, 
however, easily complicate matters by asking for the value of T for which we have 
a node occurring after exactly two turns, or three turns and so on. Following this 
logic, we would then seem to be led to the conclusion that there is something 
dramatic happening to the Kitaev chain at hnite temperature also for values of m 
which exceed 1. See Fig. |5bl It is therefore necessary to reexamine the physical 
logic behind the Uhlmann phase, and we do so in Section 


5 The interferometric geometric phase 

In Ref:s. [T^ [T3] . a different geometric phase for mixed quantum states was in¬ 
troduced in the context of interferometry. Its general dehnition is fairly involved, 
but for mixed states having non-degenerate spectra, which is the only case con¬ 
sidered here, the definition can be seen as a straightforward generalization of the 
Berry phase for pure states. 
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Consider a curve of density operators 


P(^) = (31) 

i 

such that for each k, the eigenvalues Pi{k) are non-degenerate. If the eigenstates 
develop in a parallel manner, 


{ui{k)\ui{k)) = 0 , 

we dehne the interferometric geometric phase of p{k) as 


■y{k) = argy;^/pJ(0)ft(^{tii(0)|'Ui(A^)>. 


(32) 


(33) 


As for Uhlmann’s phase, the interferometric phase does not depend on the op¬ 
erator that determines the dynamics of the system. Rather, it is a quantity 
associated with p{k) which only depends on how p{k) is embedded in the space 
of density operators. However, this is true only if the eigenstates satisfy the 


parallelism condition (32). If this condition is not met, we need to phase-shift 


the eigenstates as follows: 




(34) 


Clearly, the interferometric geometric phase reduces to the Berry phase if p{k) 
is a curve of pure states. 

The Berry phase for a cyclically evolving pure state equals half of the solid 
angle enclosed by the curve traced out by the state’s Bloch vector on the Bloch 
sphere. Using this observation, it is fairly easy to calculate the interferometric 
phase for cyclically evolving qubits, provided they do not pass the maximally 
mixed state. In fact, we believe that the parallelism condition in Eq. (32) cannot 


be extended to include all curves that pass through the maximally mixed state. 
However, this is still under investigation. 

Consider a closed curve of mixed qubits 

p{k) = pi{k)\ui{k)){ui{k)\+p 2 {k)\u 2 {k)){u 2 ik)\, (0 </c < k), (35) 


and assume that the eigenvectors \ui{k)) satisfy the parallelism condition (32) 


Let R{k) be the Bloch vector of p{k), and let f{k) = R{k)/R{k), where R{k) is 
the Euclidean length of R{k). Then f{k) is a closed curve on the Bloch sphere 
and (mi(0)|mi(/s;)) = where 9i equals half of the solid angle enclosed by f{k). 
Furthermore, ('U2(0 )|m 2(/^)) = where 62 is half of the solid angle enclosed 
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Figure 6: The nodes of the interferometric phase are situated on the radial segment 
opposite to the initial state. No build-up of the phase occurs before and after the 
passage of a node. If the curve passes the line of nodes, the phase makes a sudden 
TT-jump. 


by the curve —f{k). Accordingly, = e and the interferometric geometric 
phase of p{k) can be written as 

7 (k) = arg(pi(0)e*®^ + P2(0)e“*'^^) = arg(cos6*i -|-ii?(0) sin 6*i). (36) 


We observe, in particular, that for the planar curves dehned by ( |13[ ) and indexed 
by m, the interferometric phase is 0 if the curve has an even winding number 
with respect to the maximally mixed state, and equals vr if the winding number 
is odd, see Fig. When m = 1, however, the interferometric phase is undehned. 
Thus, it signals a closing of the band gap for the Kitaev chain in the same way 
as the Berry phase, but at a hnite temperature. 

An analysis like the one conducted for Uhlmann’s phase shows that the in¬ 
terferometric geometric phase has nodes along the radial segment in the Bloch 
ball opposite to the initial state. For the curves given by (13) this means that 
each point on the negative x-axis is a node, as shown in Fig. There is no 
phase build-up before, or after, the curve reaches a node. But at the node, the 
interferometric phase factor gets multiplied by —1. In this sense the behaviour 
of the interferometric phase closely mirrors that of the Berry phase. 


6 Discussion and conclusion 

We have been concerned with the physical interpretation, in a dehnite context, of 
two different geometric phases for mixed states—the Uhlmann and the interfer¬ 
ometric phases. Following the earlier studies El El El E], we have analyzed their 
behaviour to see if they are indicators of topological phases of matter at hnite 
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temperature. As our test case we have used the Kitaev chain, a one-dimensional 
system of non-interacting electrons which gives rise to a curve of Gibbs states in 
an equatorial plane of the Bloch ball as we move along the Brillouin zone. 

When T = 0, and the Kitaev chain is in a pure state, the Berry phase is 
a topological invariant which takes different values around m = 1 where the 
system undergoes a topological phase transition. For m < 1 the state is in a 
topologically non-trivial phase admitting localized zero-energy edge states, and 
the Berry phase factor is —1. For m > 1 there are no localized edge states 
and the Berry phase factor is -|-1. At finite temperature, on the other hand, 
the Kitaev chain is in a mixed state. Now, both of the mixed state geometric 
phase factors considered take the values ±1, but for the Uhlmann phase factor 
the behaviour is quite complex. We are not aware of a physical interpretation of 
the Uhlmann phase itself, but indeed the Uhlmann phase behaves in a way that 
may be of interest, regardless of the proposed interpretation of the states. Here 
it should be recalled that Uhlmann’s theory was designed with very different 
ends in view. It can be argued that Uhlmann’s U{N) bundle is as natural a 
construction for mixed states as is the unit sphere in Hilbert space considered as 
a U (1) bundle over the set of pure states. They both have major implications for 
the probabilistic structure of quantum mechanics HD]. For example, Uhlmann’s 
construction gives rise to a metric which is monotone under stochastic maps 
PI ES], and plays an important role in statistical inference PS]. It must also 
be mentioned that the Uhlmann phase is measurable in experiments where the 
purihed state appears as a bipartite entangled state including a system and a 
controlled ancilla PZlESlEn] — but this is a very different situation from that of 
a system in contact with a thermal bath having a huge number of uncontrollable 
degrees of freedom. 

What we found, when considering the concrete curves under study, was first 
of all an interesting, but in the context devastating, ’memory effect’. The phase 
changes along these curves occur only at certain nodes, and the locations of these 
nodes depend on how many turns around the curve have been made already. An¬ 
other observation concerns the limiting behaviour of the Uhlmann phase as the 
curve approaches the pure state boundary, which corresponds to letting the tem¬ 
perature approach zero. The Berry phase is recovered for closed curves |S], but 
for open curves this is not the case. The argument is delicate because Uhlmann’s 
bundle does not extend over the states of less than maximal rank (which play a 
special role in statistical inference). 

We now turn to the relevance of the Uhlmann phase for the condensed matter 
application we have in mind. Our first observation is that the Uhlmann phase is 
insensitive to the closure of the band gap. More precisely, if T > 0 and we let m 
approach 1 from below, we see that the Uhlmann phase changes its value before 
we reach m = 1, see Fig. When m = 1, the band gap is closed, and the curve 
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traversed by p{k) passes through the maximally mixed state. However, this does 
not affect Uhlmann’s geometric phase. In standard discussions of the topological 
origin of the edge states such behaviour of an invariant is explicitly forbidden 
(see, for instance. Ref. [301 )• Furthermore, Fig. 5b shows that the Uhlmann phase 


behaves in a very non-trivial fashion for a range of temperatures also when m > 1, 
in which case the chemical potential is in the non-topological range. We therefore 
agree with Viyuela et al. [3] when they state that the Uhlmann phase “does not 
determine the fate of the edge modes at finite temperature”. The Uhlmann phase 
was simply designed for a different purpose. 

In the spirit of Budich and Diehl [7j, we have also looked at a different geo¬ 
metric phase for mixed states, one that is sensitive to changes of the multiplicities 
in the spectra of the states. For the interferometric geometric phase, the phase 
factor is —1 for odd number of revolutions in the Brillouin zone and -1-1 for even 
numbers, no matter what the temperature is. This means that the interfero¬ 
metric phase does not detect any phase transition in temperature, but it detects 
the same phase transition as the Berry phase for pure states, i.e. at m = 1. In 
fact, it extends the Berry phase to finite temperature in a way that captures the 
periodicity of the state of the Kitaev chain. For p{k) given by Eq. (13), modeling 
pure states if T = 0 and mixed states if T > 0, induces a homomorphism between 
the fundamental groups of the Brillouin zone with end points identified and the 
equatorial plane punctured at the maximally mixed state. Both of these groups 
equal the group of integers, Z, and the homomorphism is n —)■ dn where d is the 
degree of the Bloch vector. (Thus, d = 1 if m < 1 and d = 0 if m > 1.) We 
can extend this with a second homomorphism Z —)■ Z 2 sending n to . Here 
Z 2 denotes the multiplicative group of the two elements ±1. The composite ho¬ 
momorphism n —)■ is the topological invariant given by the Berry phase for 

T = 0 and the interferometric phase for T > 0. Observe that the non-periodicity, 
due to the memory effect mentioned above, excludes the possibility of extracting 
a similar invariant from Uhlmann’s geometric phase. 

To summarize, the Uhlmann geometric phase is part of a general construction 
of great importance and import in quantum mechanics, but we see no reason for 
why it should be useful as an indicator of topological order in models resembling 
the Kitaev chain. The interferometric geometric phase on the other hand seems 
to be more appropriate as a topological order parameter at finite temperature, 
and we think it merits further investigation in the context of condensed matter. 
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